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1. Introduction 

1 According to general relativity theory, compact concentrations of energy 
(e.g., neutron stars and black holes) should warp spacetime strongly, and 
whenever such an energy concentration changes shape, it should create a 
dynamically changing spacetime warpage that propagates out through the 
Universe at the speed of light. This propagating warpage is called gravita- 
tional radiation — a name that arises from general relativity's description of 
gravity as a consequence of spacetime warpage. 

There are a number of efforts, worldwide, to detect gravitational radi- 
ation. These efforts are driven in part by the desire to "see gravitational 
waves in the flesh," but more importantly by the goal of using the waves 
as a probe of the Universe and of the nature of gravity. They should be a 
powerful probe, since they carry very detailed information about gravity 
and their sources. 

In this lecture I shall describe the prospects to study gravitational 
waves from astronomical systems that exist in our Universe today. Such 
systems are expected to radiate in the "high-frequency" (1 — 10 4 Hz) and 
"low-frequency" (10 -4 — 1Hz) gravitational-wave bands. By contrast, grav- 
itational waves from the early universe (which I shall not discuss here) 
should populate a far wider band of frequencies, ~ 10~ 18 Hz — 10 +8 Hz. 

The high-frequency band is the domain of Earth-based gravitational- 
wave detectors: resonant-mass antennas, and most especially laser inter- 
ferometers. A number of interesting compact sources fall in this band: the 
stellar collapse to a neutron star or black hole in our Galaxy and distant 



1 This paper will be published in Proceedings of IAU Symposium 165, Compact Stars 
in Binaries, edited by J. van Paradijs, E. van den Heuvel, and E. Kuulkers (Kluwer 
Academic Publishers). 
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Figure 1. The lines of force associated with the two polarizations of a gravitational 
wave. (From Ref. [1].) 



galaxies, which triggers supernovae; the rotation and vibration of neutron 
stars in our Galaxy; and the coalescence of neutron-star and stellar-mass 
black-hole binaries (M <; 1000M Q ) in distant galaxies. 

The low-frequency band is the domain of detectors flown in space (in 
Earth orbit or in interplanetary orbit). The most important of these are 
the Doppler tracking of spacecraft via microwave signals sent from Earth 
to the spacecraft and there transponded back to Earth (a technique that 
NASA has pursued since the early 1970's), and optical tracking of space- 
craft by each other (laser interferometry in space, a technique now un- 
der development for possible flight in ~ 2014). The low-frequency band 
should be populated by waves from short-period binary stars in our own 
Galaxy (main-sequence binaries, cataclysmic variables, white-dwarf bina- 
ries, neutron-star binaries, ...); by waves from white dwarfs, neutron stars, 
and small black holes spiraling into massive black holes (M ~3x 10 5 to 
3 x 1O 7 M ) in distant galaxies; and by waves from the inspiral and co- 
alescence of massive black-hole binaries (M ~ 100 to 10 8 M Q ) in distant 
galaxies. 

The body of this lecture consists of four principal sections. Section 2 
describes a network of high-frequency, ground-based laser interferometers 
that is currently under construction, and then Sec. 3 describes the high- 
frequency sources the network will search for, and the information that 
we hope to glean from their waves. Section 4 describes a system of low- 
frequency, space-based interferometers planned for launch in ~ 2014, and 
Sec. 5 describes the low-frequency sources they will search for and the 
information carried by their waves. 
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Figure 2. Schematic diagram of a laser interferometer gravitational wave detector. (From 
Ref. [1].) 

2. Ground-Based Laser Interferometers 

2.1. WAVE POLARIZATIONS, WAVEFORMS, AND HOW AN 
INTERFEROMETER WORKS 

According to general relativity theory, a gravitational wave has two lin- 
ear polarizations, conventionally called + (plus) and x (cross). Associated 
with each polarization there is a gravitational-wave field, h+ or h x , which 
oscillates in time and propagates with the speed of light. Each wave field 
produces tidal forces (gravitational stretching and squeezing forces) on any 
object or detector through which it passes. If the object is small compared 
to the waves' wavelength (as is the case for ground-based interferometers 
and resonant mass antennas), then relative to the object's center, the forces 
have the quadrupolar patterns shown in Figure 1. The names "plus" and 
"cross" are derived from the orientations of the axes that characterize the 
force patterns [2]. 

A laser interferometer gravitational wave detector ("interferometer" for 
short) consists of four masses that hang from vibration-isolated supports 
as shown in Figure 2, and the indicated optical system for monitoring the 
separations between the masses [2, 1]. Two masses are near each other, at 
the corner of an "L" , and one mass is at the end of each of the L's long 
arms. The arm lengths are nearly equal, L\ ~ L2 = L. When a gravitational 
wave, with frequencies high compared to the masses' ~ 1 Hz pendulum 
frequency, passes through the detector, it pushes the masses back and forth 
relative to each other as though they were free from their suspension wires, 
thereby changing the arm-length difference, AL = L\ — L2. That change 
is monitored by laser interferometry in such a way that the output of the 
photodiode (the interferometer's output) is directly proportional to AL(t). 
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If the gravitational waves are coming from overhead or underfoot and 
the axes of the + polarization coincide with the arms' directions, then it 
is the waves' + polarization that drives the masses, and AL(t)/L = h + (t). 
More generally, the interferometer's output is a linear combination of the 
two wave fields: 

= F + h + (t) + F x h x (t) = h(t) . (1) 

The coefficients F + and F x are of order unity and depend in a quadrupolar 
manner on the direction to the source and the orientation of the detector 
[2]. The combination h(t) of the two waveforms is called the gravitational- 
wave strain that acts on the detector; and the time evolutions of h(t), h + (t), 
and h x {t) are sometimes called waveforms. 

2.2. WAVE STRENGTHS AND INTERFEROMETER ARM LENGTHS 

The strengths of the waves from a gravitational-wave source can be esti- 
mated using the "Newtonian/quadrupole" approximation to the Einstein 
field equations. This approximation says that h ~ (G/c A )Q/r, where Q 
is the second time derivative of the source's quadrupole moment, r is the 
distance of the source from Earth, and G and c are Newton's gravitation 
constant and the speed of light. The strongest sources will be highly non- 
spherical and thus will have Q ~ ML 2 , where M is their mass and L their 
size, and correspondingly will have Q ~ 2Mv 2 ~ 4E™ Q , where v is their 
internal velocity and E£f n is the nonspherical part of their internal kinetic 
energy. This provides us with the estimate 

c z r 

i.e., h is about 4 times the gravitational potential produced at Earth by 
the mass-equivalent of the source's nonspherical, internal kinetic energy — 
made dimensionless by dividing by c 2 . Thus, in order to radiate strongly, 
the source must have a very large, nonspherical, internal kinetic energy. 

The best known way to achieve a huge internal kinetic energy is via grav- 
ity; and by energy conservation (or the virial theorem) , any gravitationally- 
induced kinetic energy must be of order the source's gravitational poten- 
tial energy. A huge potential energy, in turn, requires that the source be 
very compact, not much larger than its own gravitational radius. Thus, the 
strongest gravity- wave sources must be highly compact, dynamical concen- 
trations of large amounts of mass (e.g., colliding and coalescing black holes 
and neutron stars). 

Such sources cannot remain highly dynamical for long; their motions 
will be stopped by energy loss to gravitational waves and/or the formation 
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of an all-encompassing black hole. Thus, the strongest sources should be 
transient. Moreover, they should be very rare — so rare that to see a rea- 
sonable event rate will require reaching out through a substantial fraction 
of the Universe. Thus, just as the strongest radio waves arriving at Earth 
tend to be extragalactic, so also the strongest gravitational waves are likely 
to be extragalactic. 

For highly compact, dynamical objects that radiate in the high-frequency 
band, e.g. colliding and coalescing neutron stars and stellar-mass black 
holes, the internal, nonspherical kinetic energy -E£f n /c 2 is of order the mass 
of the Sun; and, correspondingly, Eq. (2) gives h ~ 10 -22 for such sources 
at the Hubble distance, h ~ 1CT 21 at 200 Mpc (a best-guess distance for 
several neutron-star coalescences per year; Sec. 3.1.2), and h ~ 10" 20 at the 
Virgo cluster (15 Mpc). These numbers set the scale of sensitivities that 
ground-based interferometers seek to achieve: h <J 10~ 21 to 10 -22 . 

When one examines the technology of laser interferometry, one sees good 
prospects to achieve measurement accuracies AL ~ 10~ 16 cm (1/1000 the 
diameter of the nucleus of an atom). With such an accuracy, an interfer- 
ometer must have an arm length L = AL/h ~ 1 to 10 km, in order to 
achieve the desired wave sensitivities, 10 -21 to 10~ 22 . This sets the scale of 
the interferometers that are now under construction. 

2.3. LIGO, VIRGO, AND THE INTERNATIONAL INTERFEROMETRIC 
NETWORK 

Interferometers are plagued by non-Gaussian noise, e.g. due to sudden 
strain releases in the wires that suspend the masses. This noise prevents 
a single interferometer, by itself, from detecting with confidence short- 
duration gravitational-wave bursts (though it might be possible for a single 
interferometer to search for the periodic waves from known pulsars). The 
non-Gaussian noise can be removed by cross correlating two, or preferably 
three or more, interferometers that are networked together at widely sepa- 
rated sites. 

The technology and techniques for such interferometers have been un- 
der development for nearly 25 years, and plans for km-scale interferometers 
have been developed over the past 13 years. An international network con- 
sisting of three km-scale interferometers, at three widely separated sites, is 
now in the early stages of construction. It includes two sites of the American 
LIGO Project ("Laser Interferometer Gravitational Wave Observatory") 
[1], and one site of the French/Italian VIRGO Project (named after the 
Virgo cluster of galaxies) [3]. 

LIGO will consist of two vacuum facilities with 4-kilometer-long arms, 
one in Hanford, Washington (in the northwestern United States) and the 
other in Livingston, Louisiana (in the southeastern United States). These 
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facilities are designed to house many successive generations of interferom- 
eters without the necessity of any major facilities upgrade; and after a 
planned future expansion, they will be able to house several interferometers 
at once, each with a different optical configuration optimized for a differ- 
ent type of wave (e.g., broad-band burst, or narrow-band periodic wave, 
or stochastic wave). The LIGO facilities and their first interferometers are 
being constructed by a team of about 80 physicists and engineers at Cal- 
tech and MIT, led by Barry Barish (the PI) and Gary Sanders (the Project 
Manager). Substantial contributions are also being made by scientists at 
other institutions. 

The VIRGO Project is building one vacuum facility in Pisa, Italy, with 
3-kilometer-long arms. This facility and its first interferometers are a col- 
laboration of more than a hundred physicists and engineers at the INFN 
(Frascati, Napoli, Perugia, Pisa), LAL (Orsay), LAPP (Annecy), LOA 
(Palaiseau), IPN (Lyon), ESPCI (Paris), and the University of Illinois (Ur- 
bana), under the leadership of Alain Brillet and Adalberto Giazotto. 

Both LIGO and VIRGO are scheduled for completion in the late 1990s, 
and their first gravitational-wave searches are likely to be performed in 2000 
or 2001. 

LIGO alone, with its two sites which have parallel arms, will be able 
to detect an incoming gravitational wave, measure one of its two wave- 
forms, and (from the time delay between the two sites) locate its source 
to within a ~ 1° wide annulus on the sky. LIGO and VIRGO together, 
operating as a coordinated international network, will be able to locate the 
source (via time delays plus the interferometers' beam patterns) to within 
a 2-dimensional error box with size between several tens of arcminutes and 
several degrees, depending on the source direction and on the amount of 
high-frequency structure in the waveforms. They will also be able to mon- 
itor both waveforms h+(t) and h x (t) (except for frequency components 
above about 1kHz and below about 10 Hz, where the interferometers' noise 
becomes severe). 

The accuracies of the direction measurements and the ability to monitor 
more than one waveform will be severely compromised when the source lies 
anywhere near the plane formed by the three LIGO/VIRGO interferometer 
locations. To get good all-sky coverage will require a fourth interferometer 
at a site far out of that plane; Japan and Australia would be excellent loca- 
tions, and research groups there are carrying out research and development 
on interferometric detectors, aimed at such a possibility. A 300-meter pro- 
totype interferometer called TAMA is under construction in Tokyo, and a 
400-meter prototype called AIGO400 has been proposed for construction 
north of Perth. 

Two other groups are major players in this field, one in Britain led by 
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James Hough, the other in Germany, led by Karsten Danzmann. These 
groups each have two decades of experience with prototype interferometers 
(comparable experience to the LIGO team and far more than anyone else) 
and great expertise. Frustrated by inadequate financing for a kilometer- 
scale interferometer, they are constructing, instead, a 600 meter system 
called GEO600 near Hannover, Germany. Their goal is to develop, from 
the outset, an interferometer with the sort of advanced design that LIGO 
and VIRGO will attempt only as a "second-generation" instrument, and 
thereby achieve sufficient sensitivity to be full partners in the international 
network's first gravitational-wave searches; they then would offer a variant 
of their interferometer as a candidate for second-generation operation in 
the much longer arms of LIGO and/or VIRGO. It is a seemingly audacious 
plan, but with their extensive experience and expertise, the British/German 
collaboration might pull it off successfully. 

3. High-Frequency Gravitational- Wave Sources 

3.1. COALESCING COMPACT BINARIES 

The best understood of all gravitational-wave sources are coalescing, com- 
pact binaries composed of neutron stars (NS) and black holes (BH). These 
NS/NS, NS/BH, and BH/BH binaries may well become the "bread and 
butter" of the LIGO/VIRGO diet. 

The Hulse- Taylor [4, 5] binary pulsar, PSR 1913+16, is an example 
of a NS/NS binary whose waves could be measured by LIGO/VIRGO, 
if we were to wait long enough. At present PSR1913+16 has an orbital 
frequency of about 1/(8 hours) and emits its waves predominantly at twice 
this frequency, roughly 10~ 4 Hz, which is in the low-frequency band — far 
too low to be detected by LIGO/VIRGO. However, as a result of their 
loss of orbital energy to gravitational waves, the PSR1913+16 NS's are 
gradually spiraling inward. If we wait roughly 10 8 years, this inspiral will 
bring the waves into the LIGO/VIRGO high-frequency band. As the NS's 
continue their inspiral, the waves will then sweep upward in frequency, over 
a time of about 15 minutes, from 10 Hz to ~ 10 3 Hz, at which point the 
NS's will collide and coalesce. It is this last 15 minutes of inspiral, with 
~ 16, 000 cycles of waveform oscillation, and the final coalescence, that 
LIGO/VIRGO seeks to monitor. 

3.1.1. Wave Strengths Compared to LIGO Sensitivities 
Figure 3 compares the projected sensitivities of interferometers in LIGO [1] 
with the wave strengths from the last few minutes of inspiral of BH/BH, 
NS/BH, and NS/NS binaries at various distances from Earth. The two solid 
curves at the bottoms of the stippled regions (labeled h ims ) are the rms noise 
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Figure 3. LIGO's projected broad-band noise h Ilns and sensitivity to bursts hsB [1] 
compared with the strengths of the waves from the last few minutes of inspiral of compact 
binaries. The signal to noise ratios are V2 higher than in Ref. [1] because of a factor 2 
error in Eq. (29) of Ref. [2]. 



levels for broad-band waves that have optimal direction and polarization. 
The tops of the stippled regions (labeled hsB for "sensitivity to bursts") 
are the sensitivities for highly confident detection of randomly polarized, 
broad-band waves from random directions (i.e., the sensitivities for high 
confidence that any such observed signal is not a false alarm due to Gaussian 
noise). The upper stippled region and its bounding curves are the expected 
performances of the first interferometers in LIGO; the lower stippled region 
and curves are performances of more advanced LIGO interferometers. 

As the NS's and/or BH's spiral inward, their waves sweep upward in 
frequency (left to right in the diagram). The dashed lines show their "char- 
acteristic" signal strength h c (approximately the amplitude h of the waves' 
oscillations multiplied by the square root of the number of cycles spent near 
a given frequency, \/n); the signal-to-noise ratio is this h c divided by the 
detector's \/5h vms , S/N = h c /(\^5h ims ), where the \/5 converts h rms from 
"optimal direction and polarization" to "random direction and polariza- 
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tion") [1, 2]. The arrows along the NS/NS inspiral track indicate the time 
until final coalescence and the separation between the NS centers of mass. 
Each NS is assumed to have a mass of 1.4 suns and a radius ~ 10 km; each 
BH, 10 suns and ~ 20 km. 

Notice that the signal strengths in Fig. 3 are in good accord with rough 
estimates based on Eq. (2); at the endpoint (right end) of each inspiral, 
the number of cycles n spent near that frequency is of order unity, so the 
quantity plotted, h c ~ h^/n, is about equal to h — and at distance 200 Mpc 
is about 10~ 21 , as Eq. (2) predicts. 

3.1.2. Coalescence Rates 

Such final coalescences are few and far between in our own galaxy: about 
one every 100,000 years, according to 1991 estimates by Phinney [6] and 
by Narayan, Piran, and Shemi [7], based on the statistics of binary pulsar 
searches in our galaxy which found three that will coalesce in less than 10 10 
years. Extrapolating out through the universe on the basis of the density 
of production of blue light (the color produced predominantly by massive 
stars), Phinney [6] and Narayan et. al. [7] infer that to see several NS/NS 
coalescences per year, LIGO/VIRGO will have to look out to a distance of 
about 200 Mpc (give or take a factor ~ 2); cf. the "NS/NS inspiral, 200 
Mpc" line in Fig. 3. Since these estimates were made, the binary pulsar 
searches have been extended through a significantly larger volume of the 
galaxy than before, and no new ones with coalescence times ^ 10 10 years 
have been found; as a result, the binary-pulsar-search-based best estimate 
of the coalescence rate should be revised downward [8], perhaps to as little 
as one every million years in our galaxy, corresponding to a distance 400 
Mpc for several per year [8]. 

A rate of one every million years in our galaxy is ~ 1000 times smaller 
than the birth rate of the NS/NS binaries' progenitors: massive, compact, 
main-sequence binaries [6, 7]. Therefore, either 99.9 per cent of progeni- 
tors fail to make it to the NS/NS state (e.g., because of binary disruption 
during a supernova or forming TZO's), or else they do make it, but they 
wind up as a class of NS/NS binaries that has not yet been discovered in 
any of the pulsar searches. Several experts on binary evolution have argued 
for the latter [9, 10, 11, 12]: most NS/NS binaries, they suggest, may form 
with such short orbital periods that their lifetimes to coalescence are signif- 
icantly shorter than normal pulsar lifetimes (~ 10 years); and with such 
short lifetimes, they have been missed in pulsar searches. By modeling the 
evolution of the galaxy's binary star population, the binary experts arrive 
at best estimates as high as 3 x 10 -4 coalescences per year in our galaxy, 
corresponding to several per year out to 60 Mpc distance [9]. Phinney [6] 
describes other plausible populations of NS/NS binaries that could increase 



10 



KIP S. THORNE 



the event rate, and he argues for "ultraconservative" lower and upper limits 
of 23 Mpc and lOOOMpc for how far one must look to see several coalescence 
per year. 

By comparing these rate estimates with the signal strengths in Fig. 3, 
we see that (i) the first interferometers in LIGO/VIRGO have a moderate 
but not high probability of seeing NS/NS coalescences; (ii) advanced inter- 
ferometers are almost certain of seeing them (the requirement that this be 
so was one factor that forced the LIGO/VIRGO arm lengths to be so long, 
several kilometers); and (iii) they are most likely to be discovered roughly 
half-way between the first and advanced interferometers — which means by 
an improved variant of the first interferometers several years after LIGO 
operations begin. 

We have no good observational handle on the coalescence rate of NS/BH 
or BH/BH binaries. However, theory suggests that their progenitors might 
not disrupt during the stellar collapses that produce the NS's and BH's, so 
their coalescence rate could be about the same as the birth rate for their 
progenitors: ~ 1/100, 000 years in our galaxy. This suggests that within 200 
Mpc distance there might be several NS/BH or BH/BH coalescences per 
year. [6, 7, 9, 12]. This estimate should be regarded as a plausible upper 
limit on the event rate and lower limit on the distance to look [6, 7]. 

If this estimate is correct, then NS/BH and BH/BH binaries will be seen 
before NS/NS, and might be seen by the first LIGO/VIRGO interferometers 
or soon thereafter; cf. Fig. 3. However, this estimate is far less certain than 
the (rather uncertain) NS/NS estimates! 

Once coalescence waves have been discovered, each further improvement 
of sensitivity by a factor 2 will increase the event rate by 2 3 ~ 10. Assuming 
a rate of several NS/NS per year at 200 Mpc, the advanced interferometers 
of Fig. 3 should see ~ 100 per year. 

3.1.3. Inspiral Waveforms and the Information They Can Bring 

Neutron stars and black holes have such intense self gravity that it is ex- 
ceedingly difficult to deform them. Correspondingly, as they spiral inward 
in a compact binary, they do not gravitationally deform each other signifi- 
cantly until several orbits before their final coalescence [13, 14]. This means 
that the inspiral waveforms are determined to high accuracy by only a few, 
clean parameters: the masses and spin angular momenta of the bodies, 
and the initial orbital elements (i.e. the elements when the waves enter the 
LIGO/VIRGO band). 

Though tidal deformations are negligible during inspiral, relativistic ef- 
fects can be very important. If, for the moment, we ignore the relativistic 
effects — i.e., if we approximate gravity as Newtonian and the wave gener- 
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Figure 4- Waveforms from the inspiral of a compact binary, computed using Newtonian 
gravity for the orbital evolution and the quadrupole moment approximation for the wave 
generation. (From Ref. [1].) 



ation as due to the binary's oscillating quadrupole moment [2], then the 
shapes of the inspiral waveforms h+(t) and h x (t) are as shown in Fig. 4. 

The left-hand graph in Fig. 4 shows the waveform increasing in ampli- 
tude and sweeping upward in frequency (i.e., undergoing a "chirp") as the 
binary's bodies spiral closer and closer together. The ratio of the amplitudes 
of the two polarizations is determined by the inclination t of the orbit to our 
line of sight (lower right in Fig. 4). The shapes of the individual waves, i.e. 
the waves' harmonic content, are determined by the orbital eccentricity (up- 
per right). (Binaries produced by normal stellar evolution should be highly 
circular due to past radiation reaction forces, but compact binaries that 
form by capture events, in dense star clusters that might reside in galactic 
nuclei [15], could be quite eccentric.) If, for simplicity, the orbit is circular, 
then the rate at which the frequency sweeps or "chirps", df/dt [or equiva- 
lently the number of cycles spent near a given frequency, n = f 2 (df /dt)~ l \ 
is determined solely, in the Newtonian/quadrupole approximation, by the 
binary's so-called chirp mass, M c = (M 1 M 2 ) 3/5 /(M 1 + M2) 1 / 5 (where M x 
and M2 are the two bodies' masses). The amplitudes of the two waveforms 
are determined by the chirp mass, the distance to the source, and the orbital 
inclination. Thus (in the Newtonian/quadrupole approximation), by mea- 
suring the two amplitudes, the frequency sweep, and the harmonic content 
of the inspiral waves, one can determine as direct, resulting observables, 
the source's distance, chirp mass, inclination, and eccentricity [16, 17]. 

As in binary pulsar observations [5], so also here, relativistic effects 
add further information: they influence the rate of frequency sweep and 
produce waveform modulations in ways that depend on the binary's di- 
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Figure 5. Modulational envelope for the waveform from a IMq nonspinning NS spiraling 
into a IOMq , rapidly spinning Kerr black hole (spin parameter a = 1). The orbital angular 
momentum L is inclined by a — 11.3 degrees to the hole's spin angular momentum S, 
and the two precess around J = L + S, whose direction remains fixed in space as L = |L| 
shrinks and S — |S| = Mbhcl remains constant. The precession modulates the waves by 
an amount that depends on (i) the direction to Earth (here along the initial L x S, i.e. out 
of the paper) and (ii) the orientation of the detector's arms (here parallel to the figure's 
initial L and to Lx (direction to Earth) for h+, and rotated 45 degrees for hx)- The 
figure shows the waveforms' modulational envelopes (in arbitrary units, the same for h+ 
and h x ), parametrized by the wave frequency / and the number of cycles of oscillation 
between the indicated /'s. The total number of precessions from / to coalescence is 
7V prcc ~ (5/64n)(Ma/Li)(nMf)~ 2/3 ~ 20(//10Hz)" 2/3 . (From [18, 19].) 

mensionless ratio r\ = n/M of reduced mass \i = M1M2/ (M\ + M2) to 
total mass M = M\ + M2 [20] and on the spins of the binary's two bodies 
[21]; These relativistic effects are reviewed and discussed at length in Refs. 
[18, 22]. Two deserve special mention: (i) As the waves emerge from the 
binary, some of them get backscattered one or more times off the binary's 
spacetime curvature, producing wave tails. These tails act back on the bi- 
nary, modifying its inspiral rate in a measurable way. (ii) If the orbital 
plane is inclined to one or both of the binary's spins, then the spins drag 
inertial frames in the binary's vicinity (the "Lense-Thirring effect"), this 
frame dragging causes the orbit to precess, and the precession modulates 
the waveforms [18, 19, 23]. Figure 5 shows the resulting modulation for a 
IMq NS spiraling into a rapidly spinning, IOMq BH. 

Remarkably, the relativistic corrections to the frequency sweep will be 
measurable with very high accuracy, even though they are typically <i 10 
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per cent of the Newtonian contribution, and even though the typical signal 
to noise ratio will be only ~ 9 even after optimal signal processing. The 
reason is as follows [24, 25, 18]: 

The frequency sweep will be monitored by the method of "matched fil- 
ters"; in other words, the incoming, noisy signal will be cross correlated 
with theoretical templates. If the signal and the templates gradually get 
out of phase with each other by more than ~ 1/10 cycle as the waves sweep 
through the LIGO/VIRGO band, their cross correlation will be significantly 
reduced. Since the total number of cycles spent in the LIGO/VIRGO band 
will be ~ 16,000 for a NS/NS binary, ~ 3500 for NS/BH, and ~ 600 for 
BH/BH, this means that LIGO/VIRGO should be able to measure the 
frequency sweep to a fractional precision <^ 10~ 4 , compared to which the 
relativistic effects are very large. (This is essentially the same method as 
Joseph Taylor and colleagues use for high-accuracy radio-wave measure- 
ments of relativistic effects in binary pulsars [5].) 

Preliminary analyses, using the theory of optimal signal processing, 
predict the following typical accuracies for LIGO/VIRGO measurements 
based solely on the frequency sweep (i.e., ignoring modulational informa- 
tion) [26, 24, 25, 27, 18]: (i) The chirp mass M c will typically be measured, 
from the Newtonian part of the frequency sweep, to ~ 0.04% for a NS/NS 
binary and ~ 0.3% for a system containing at least one BH. (ii) If we are 
confident (e.g., on a statistical basis from measurements of many previous 
binaries) that the spins are a few percent or less of the maximum physi- 
cally allowed, then the reduced mass // will be measured to ~ 1% for NS/NS 
and NS/BH binaries, and ~ 3% for BH/BH binaries. (Here and below NS 
means a ~ 1.4M neutron star and BH means a ~ 1OM black hole.) (iii) 
Because the frequency dependences of the (relativistic) fx effects and spin 
effects are not sufficiently different to give a clean separation between /j, 
and the spins, if we have no prior knowledge of the spins, then the spin/^ 
correlation will worsen the typical accuracy of fx by a large factor, to ~ 30% 
for NS/NS, ~ 50% for NS/BH, and a factor ~ 2 for BH/BH [26, 24]. These 
worsened accuracies might be improved somewhat by waveform modula- 
tions caused by the spin-induced precession of the orbit [19, 23], and even 
without modulational information, a certain combination of fi and the spins 
will be determined to a few per cent. Much additional theoretical work is 
needed to firm up the measurement accuracies. 

To take full advantage of all the information in the inspiral wave- 
forms will require theoretical templates that are accurate, for given masses 
and spins, to a fraction of a cycle during the entire sweep through the 
LIGO/VIRGO band. Such templates are being computed by an interna- 
tional consortium of relativity theorists (Blanchet and Damour in France, 
Iyer in India, Will and Wiseman in the U.S., and others) [22, 28], using 
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post-Newtonian expansions of the Einstein field equations. This enterprise 
is rather like computing the Lamb shift to high order in powers of the fine 
structure constant, for comparison with experiment. The terms of leading 
order in the mass ratio rj = fi/M are being checked by a Japanese- American 
consortium (Poisson, Nakamura, Sasaki, Tagoshi, Tanaka) using the Teukol- 
sky formalism for weak perturbations of black holes [29, 30]. These small-?? 
calculations have been carried to very high post-Newtonian order for circu- 
lar orbits and no spins [31, 32], and from those results Cutler and Flanagan 
[33] have estimated the order to which the full, finite-?? computations must 
be carried in order that systematic errors in the theoretical templates will 
not significantly impact the information extracted from the LIGO/VIRGO 
observational data. The answer appears daunting: radiation-reaction ef- 
fects must be computed to three full post-Newtonian orders [six orders in 
v/c =(orbital velocity)/(speed of light)] beyond the leading-order radiation 
reaction, which itself is 5 orders in v/c beyond the Newtonian theory of 
gravity. 

It is only about ten years since controversies over the leading-order radi- 
ation reaction [34] were resolved by a combination of theoretical techniques 
and binary pulsar observations. Nobody dreamed then that LIGO/VIRGO 
observations will require pushing post-Newtonian computations onward 
from 0[(w/c) 5 ] to 0[(v / c) 11 ]. This requirement epitomizes a major change 
in the field of relativity research: At last, 80 years after Einstein formulated 
general relativity, experiment has become a major driver for theoretical 
analyses. 

Remarkably, the goal of 0[(v/c) n ] is achievable. The most difficult part 
of the computation, the radiation reaction, has been evaluated to 0[{v/cf\ 
beyond Newton by the French/Indian/ American consortium [28] and as of 
this writing, rumors have it that 0[(v/c) 10 ] is coming under control. 

These high-accuracy waveforms are needed only for extracting informa- 
tion from the inspiral waves, after the waves have been discovered; they 
are not needed for the discovery itself. The discovery is best achieved us- 
ing a different family of theoretical waveform templates, one that covers the 
space of potential waveforms in a manner that minimizes computation time 
instead of a manner that ties quantitatively into general relativity theory 
[18]. Such templates are in the early stage of development [35, 36, 37]. 

LIGO/VIRGO observations of compact binary inspiral have the poten- 
tial to bring us far more information than just binary masses and spins: 
— They can be used for high-precision tests of general relativity. In scalar- 
tensor theories (some of which are highly attractive alternatives to 
general relativity [38]), radiation reaction due to emission of scalar 
waves places a unique signature on those waves that LIGO/VIRGO 
would detect — a signature that can be searched for with high precision 
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[39]. 

— They can be used to measure the Hubble constant, deceleration param- 
eter, and cosmological constant [16, 17, 40, 41]. The keys to such mea- 
surements are that (i) advanced interferometers in LIGO/ VIRGO will 
be able to see NS/NS out to cosmological redshifts z ~ 0.3, and NS/BH 
out to z ~ 2. (ii) The direct observables that can be extracted from 
the observed waves include the source's luminosity distance tl (mea- 
sured to accuracy ~ 10 per cent in a large fraction of cases), and its 
direction on the sky (to accuracy ~ 1 square degree) — accuracies good 
enough that only one or a few electromagnetically-observed clusters of 
galaxies should fall within the 3-dimensional gravitational error boxes, 
thereby giving promise to joint gravitational/electromagnetic statisti- 
cal studies, (iii) Another direct gravitational observable is (1 + z)M 
where z is redshift and M is any mass in the system (measured to the 
accuracies quoted above). Since the masses of NS's in binaries seem to 
cluster around 1.4M , measurements of (l + z)M can provide a handle 
on the redshift, even in the absence of electromagnetic aid. 

— For a NS or small BH spiraling into a massive ~ 50 to 500-M© BH, the 
inspiral waves will carry a "map" of the spacetime geometry around 
the big hole — a map that can be used, e.g., to test the theorem that 
"a black hole has no hair" [42]; cf. Sec. 5.3 below. 

3.1.4. Coalescence Waveforms and their Information 

The waves from the binary's final coalescence can bring us new types of 
information. 

BH/BH Coalescence: In the case of a BH/BH binary, the coalescence 
will excite large-amplitude, highly nonlinear vibrations of spacetime curva- 
ture near the coalescing black-hole horizons — a phenomenon of which we 
have very little theoretical understanding today. Especially fascinating will 
be the case of two spinning black holes whose spins are not aligned with 
each other or with the orbital angular momentum. Each of the three angular 
momentum vectors (two spins, one orbital) will drag space in its vicinity 
into a tornado-like swirling motion — the general relativistic "dragging of 
inertial frames," so the binary is rather like two tornados with orientations 
skewed to each other, embedded inside a third, larger tornado with a third 
orientation. The dynamical evolution of such a complex configuration of 
spacetime warpage (as revealed by its emitted waves) may well bring us 
surprising new insights into relativistic gravity. Moreover, if the sum of the 
BH masses is fairly large, ~ 40 to 2OOM , then the waves should come off 
in a frequency range / ~ 40 to 200 Hz where the LIGO/VIRGO broad- 
band interferometers have their best sensitivity and can best extract the 
information the waves carry. 
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To get full value out of such wave observations will require having the- 
oretical computations with which to compare them [43]. There is no hope 
to perform such computations analytically; they can only be done as su- 
percomputer simulations. The development of such simulations is a major 
effort within the world's relativity community. 

NS/NS Coalescence: The final coalescence of NS/NS binaries should 
produce waves that are sensitive to the equation of state of nuclear matter, 
so such coalescences have the potential to teach us about the nuclear equa- 
tion of state [18]. In essence, we will be studying nuclear physics via the 
collisions of atomic nuclei that have nucleon numbers A ~ 10 57 — somewhat 
larger than physicists are normally accustomed to. The accelerator used to 
drive these nuclei up to the speed of light is the binary's self gravity, and the 
radiation by which the details of the collisions are probed is gravitational. 

A number of research groups [44, 45, 46, 47, 13, 14, 48, 49] are engaged 
in numerical astrophysics simulations of NS/NS coalescence, with the goal 
not only to predict the emitted gravitational waveforms and their depen- 
dence on equation of state, but also (more immediately) to learn whether 
such coalescences might power the 7-ray bursts that have been a major 
astronomical puzzle since their discovery in the early 1970s. If advanced 
LIGO interferometers were now in operation, they could report definitively 
whether or not the 7-bursts are produced by NS/NS binaries; and if the 
answer were yes, then the combination of 7-burst data and gravitational- 
wave data could bring valuable information that neither could bring by 
itself. For example, we could determine when, to within a few msec, the 
7-burst is emitted relative to the moment the NS's first begin to touch; 
and by comparing the 7 and gravitational times of arrival, we might test 
whether gravitational waves propagate with the speed of light to a frac- 
tional precision of ~ 0.01sec/3 x 10 9 lyr ~ 1(T 19 . 

Unfortunately, the final NS/NS coalescence will emit its gravitational 
waves in the kHz frequency band (800Hz <^ / <^ 2500Hz) where pho- 
ton shot noise will prevent them from being studied by the standard, 
"workhorse," broad-band interferometers of Fig. 3. However, a specially 
configured ("dual-recycled") interferometer invented by Brian Meers [50], 
which could have enhanced sensitivity in the kHz region at the price of 
reduced sensitivity elsewhere, may be able to measure the waves and ex- 
tract their equation of state information, as might massive, spherical bar 
detectors [18, 51]. Such measurements will be very difficult and are likely 
only when the LIGO/VIRGO network has reached a mature stage. 

NS/BH Coalescence: A NS spiraling into a BH of mass M £ 1OM 
should be swallowed more or less whole. However, if the BH is less massive 
than roughly 1OM , and especially if it is rapidly rotating, then the NS will 
tidally disrupt before being swallowed. Little is known about the disruption 
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and accompanying waveforms. To model them with any reliability will likely 
require full numerical relativity, since the circumferences of the BH and 
NS will be comparable and their physical separation at the moment of 
disruption will be of order their separation. As with NS/NS, the coalescence 
waves should carry equation of state information and will come out in the 
kHz band, where their detection will require advanced, specialty detectors. 

Christodoulou Memory: As the coalescence waves depart from their 
source, their energy creates (via the nonlinearity of Einstein's field equa- 
tions) a secondary wave called the "Christodoulou memory" [52, 53, 54]. 
Whereas the primary waves may have frequencies in the kHz band, the 
memory builds up on the timescale of the primary energy emission pro- 
file, which is likely to be of order 0.01 sec, corresponding to a memory 
frequency in the optimal band for the LIGO/VIRGO workhorse interfer- 
ometers, ~ 100Hz. Unfortunately, the memory is so weak that only very 
advanced interferometers have much chance of detecting and studying it — 
and then, perhaps only for BH/BH coalescences and not for NS/NS or 
NS/BH [55]. 

3.2. STELLAR CORE COLLAPSE AND SUPERNOVAE 

Several features of the stellar core collapse, which triggers supernovae, can 
produce significant gravitational radiation in the high-frequency band. We 
shall consider these features in turn, the most weakly radiating first. 

3.2.1. Boiling of the Newborn Neutron Star 

Even if the collapse is spherical, so it cannot radiate any gravitational waves 
at all, it should produce a convectively unstable neutron star that "boils" 
vigorously (and nonspherically) for the first ~ 0.1 second of its life [56]. The 
boiling dredges up high-temperature nuclear matter (T ~ 10 12 K) from the 
neutron star's central regions, bringing it to the surface (to the "neutrino- 
sphere" ) , where it cools by neutrino emission before being swept back down- 
ward and reheated. Burrows estimates [57, 58] that the boiling should gen- 
erate n ~ 10 cycles of gravitational waves with frequency / ~ 100Hz and 
amplitude h ~ 3 x 10~ 22 (30kpc/r) (where r is the distance to the source), 
corresponding to a characteristic amplitude h c ~ h^fn ~ 10~ 21 (30kpc/r); 
cf. Fig. 6. LIGO/VIRGO will be able to detect such waves only in the lo- 
cal group of galaxies, where the supernova rate is probably no larger than 
~ 1 each 10 years. However, neutrino detectors have a similar range, and 
there could be a high scientific payoff from correlated observations of the 
gravitational waves emitted by the boiling's mass motions and neutrinos 
emitted from the boiling neutrino-sphere. 
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Figure 6. Characteristic amplitudes of the gravitational waves from various processes 
accompanying stellar core collapse and supernovae, compared with projected sensitivities 
of LIGO's interferometers. 



3.2.2. Axisymmetric Collapse, Bounce, and Oscillations 

Rotation will centrifugally flatten the collapsing core, enabling it to radi- 
ate as it implodes. If the core's angular momentum is small enough that 
centrifugal forces do not halt or strongly slow the collapse before it reaches 
nuclear densities, then the core's collapse, bounce, and subsequent oscil- 
lations are likely to be axially symmetric. Numerical simulations [59, 60] 
show that in this case the waves from collapse, bounce, and oscillation will 
be quite weak: the total energy radiated as gravitational waves is not likely 
to exceed ~ 1CP 7 solar masses (about 1 part in a million of the collapse 
energy) and might often be much less than this; and correspondingly, the 
waves' characteristic amplitude will be h c <^ 3 x lCT 21 (30kpc/r). These 
collapse-and-bounce waves will come off at frequencies ~ 200 Hz to ~ 1000 
Hz, and will precede the boiling waves by a fraction of a second. Like the 
boiling waves, they probably cannot be seen by LIGO/VIRGO beyond the 
local group of galaxies and thus will be a very rare occurrence. 
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3.2.3. Rotation- Induced Bars and Break-Up 

If the core's rotation is large enough to strongly flatten the core before or 
as it reaches nuclear density, then a dynamical and/or secular instability is 
likely to break the core's axisymmetry. The core will be transformed into 
a bar-like configuration that spins end-over-end like an American football, 
and that might even break up into two or more massive pieces. In this case, 
the radiation from the spinning bar or orbiting pieces could be almost as 
strong as that from a coalescing neutron-star binary, and thus could be seen 
by the LIGO/ VIRGO first interferometers out to the distance of the Virgo 
cluster (where the supernova rate is several per year) and by advanced 
interferometers out to several hundred Mpc (supernova rate ~ 10 4 per 
year); cf. Fig. 6. It is far from clear what fraction of collapsing cores will 
have enough angular momentum to break their axisymmetry, and what 
fraction of those will actually radiate at this high rate; but even if only 
~ 1/1000 or 1/10 4 do so, this could ultimately be a very interesting source 
for LIGO/VIRGO. 

Several specific scenarios for such non-axisymmetry have been identified: 
Centrifugal hangup at ~ 100km radius: If the pre-collapse core is 
rapidly spinning (e.g., if it is a white dwarf that has been spun up by accre- 
tion from a companion), then the collapse may produce a highly flattened, 
centrifugally supported disk with most of its mass at radii R ~ 100km, 
which then (via instability) may transform itself into a bar or may bifurcate. 
The bar or bifurcated lumps will radiate gravitational waves at twice their 
rotation frequency, / ~ 100Hz — the optimal frequency for LIGO/VIRGO 
interferometers. To shrink on down to ~ 10km size, this configuration must 
shed most of its angular momentum. // a substantial fraction of the an- 
gular momentum goes into gravitational waves, then independently of the 
strength of the bar, the waves will be nearly as strong as those from a 
coalescing binary. The reason is this: The waves' amplitude h is propor- 
tional to the bar's ellipticity e, the number of cycles n of wave emission is 
proportional to 1/e 2 , and the characteristic amplitude h c = hy/n is thus 
independent of the ellipticity and is about the same whether the configu- 
ration is a bar or is two lumps [17]. The resulting waves will thus have h c 
roughly half as large, at / ~ 100Hz, as those from a NS/NS binary (half as 
large because each lump might be half as massive as a NS), and they will 
chirp upward in frequency in a manner similar to those from a binary. 

It is rather likely, however, that most of the excess angular momentum 
does not go into gravitational waves, but instead goes largely into hydro- 
dynamic waves as the bar or lumps, acting like a propeller, stir up the 
surrounding stellar mantle. In this case, the radiation will be correspond- 
ingly weaker. 

Centrifugal hangup at R ~ 20km: Lai and Shapiro [61] have ex- 
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plored the case of centrifugal hangup at radii not much lager than the final 
neutron star, say R ~ 20km. Using compressible ellipsoidal models, they 
have deduced that, after a brief period of dynamical bar-mode instability 
with wave emission at / ~ 1000Hz (explored by Houser, Centrella, and 
Smith [62]), the star switches to a secular instability in which the bar's 
angular velocity gradually slows while the material of which it is made re- 
tains its high rotation speed and circulates through the slowing bar. The 
slowing bar emits waves that sweep downward in frequency through the 
LIGO/VIRGO optimal band / ~ 100Hz, toward ~ 10Hz. The characteristic 
amplitude (Fig. 6) is only modestly smaller than for the upward-sweeping 
waves from hangup at R ~ 100km, and thus such waves should be de- 
tectable near the Virgo Cluster by the first LIGO/VIRGO interferometers, 
and at distances of a few lOOMpc by advanced interferometers. 

Successive fragmentations of an accreting, newborn neutron 
star: Bonnell and Pringle [63] have focused on the evolution of the rapidly 
spinning, newborn neutron star as it quickly accretes more and more mass 
from the pre-supernova star's inner mantle. If the accreting material carries 
high angular momentum, it may trigger a renewed bar formation, lump 
formation, wave emission, and coalescence, followed by more accretion, bar 
and lump formation, wave emission, and coalescence. Bonnell and Pringle 
speculate that hydrodynamics, not wave emission, will drive this evolution, 
but that the total energy going into gravitational waves might be as large 
as ~ 1O~ 3 M . This corresponds to h c ~ 10~ 21 (10Mpc/r). 

3.3. SPINNING NEUTRON STARS; PULSARS 

As the neutron star settles down into its final state, its crust begins to 
solidify (crystalize) . The solid crust will assume nearly the oblate axisym- 
metric shape that centrifugal forces are trying to maintain, with poloidal 
ellipticity e p oc (angular velocity of rotation) 2 . However, the principal axis 
of the star's moment of inertia tensor may deviate from its spin axis by 
some small "wobble angle" 9 W , and the star may deviate slightly from ax- 
isymmetry about its principal axis; i.e., it may have a slight ellipticity e e 
in its equatorial plane. 

As this slightly imperfect crust spins, it will radiate gravitational waves 
[64]: e e radiates at twice the rotation frequency, / = 2/ rot with h oc e e , 
and the wobble angle couples to e p to produce waves at / = / rot + / pre c 
(the precessional sideband of the rotation frequency) with amplitude h oc 
6 w e p . For typical neutron-star masses and moments of inertia, the wave 
amplitudes are 

fe~6x 10~ 25 (Jg-) 2 (^) ( ee Z e f P ) • (3) 
V500H Z y \ r J \ lO- 6 J v ' 
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The neutron star gradually spins down, due in part to gravitational- 
wave emission but perhaps more strongly due to electromagnetic torques 
associated with its spinning magnetic field and pulsar emission. This spin- 
down reduces the strength of centrifugal forces, and thereby causes the 
star's poloidal ellipticity e p to decrease, with an accompanying breakage 
and resolidification of its crust's crystal structure (a "starquake") [65]. In 
each starquake, 9 W , e e , and e p will all change suddenly, thereby changing 
the amplitudes of the star's two gravitational "spectral lines" / = 2/ rot 
and / = / rot + /prec- After each quake, there should be a healing period in 
which the star's fluid core and solid crust, now rotating at different speeds, 
gradually regain synchronism. By monitoring the amplitudes, frequencies, 
and phases of the two gravitational-wave spectral lines, and by comparing 
with timing of the electromagnetic pulsar emission, one might learn much 
about the physics of the neutron-star interior. 

How large will the quantities e e and 9 w e p be? Rough estimates of the 
crustal shear moduli and breaking strengths suggest an upper limit in the 
range e max ~ 10~ 4 to 10~ 6 , and it might be that typical values are far below 
this. We are extremely ignorant, and correspondingly there is much to be 
learned from searches for gravitational waves from spinning neutron stars. 

One can estimate the sensitivity of LIGO /VIRGO (or any other broad- 
band detector) to the periodic waves from such a source by multiplying the 
waves' amplitude h by the square root of the number of cycles over which 
one might integrate to find the signal, n = ff where f is the integration 
time. The resulting effective signal strength, h^fn, is larger than h by 

/ f \ 1/2 / f \ 1/2 

^=y7f = io 5 (_ I — ) ( — - — ) . (4) 

v v J VIOOOHzy V4month S y v ' 

This hi/n should be compared (i) to the detector's rms broad-band noise 
level for sources in a random direction, v5ft, rms , to deduce a signal-to- noise 
ratio, or (ii) to hgB to deduce a sensitivity for high-confidence detection 
when one does not know the waves' frequency in advance [2]. Such a com- 
parison suggests that the first interferometers in LIGO/VIRGO might pos- 
sibly see waves from nearby spinning neutron stars, but the odds of success 
are very unclear. 

The deepest searches for these nearly periodic waves will be performed 
by narrow-band detectors, whose sensitivities are enhanced near some cho- 
sen frequency at the price of sensitivity loss elsewhere — e.g., "dual recy- 
cled" interferometers or resonant bars. With "advanced-detector technol- 
ogy," dual-recycled interferometers might be able to detect with confidence 
all spinning neutron stars that have [2] 

(£eOr9 „ £p) ^3xl0-.»(^) 2 (^) 2 . (5) 
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There may well be a large number of such neutron stars in our galaxy; but 
it is also conceivable that there are none. We are extremely ignorant. 

Some cause for optimism arises from several physical mechanisms that 
might generate radiating ellipticities large compared to 3 x 10~ 10 : 

— It may be that, inside the superconducting cores of many neutron stars, 
there are trapped magnetic fields with mean strength B core ~ 10 13 G or 
even 10 15 G. Because such a field is actually concentrated in flux tubes 
with B = -Bcrit ~ 6 x 10 14 G surrounded by field-free superconductor, 
its mean pressure is ps = B core B cr n/8ir. This pressure could produce 
a radiating ellipticity e e ~ 6 w e p ~ Pb/p ~ 10~ 8 _B CO re/10 13 G (where p 
is the core's material pressure). 

— Accretion onto a spinning neutron star can drive precession (keeping 9 W 
substantially nonzero), and thereby might produce measurably strong 
waves [66]. 

— If a neutron star is born rotating very rapidly, then it may experience a 
gravitational-radiation-reaction-driven instability. In this "CFS" (Chan- 
drasekhar, [67] Friedman, Schutz [68]) instability, density waves prop- 
agate around the star in the opposite direction to its rotation, but are 
dragged forward by the rotation. These density waves produce gravi- 
tational waves that carry positive energy as seen by observers far from 
the star, but negative energy from the star's viewpoint; and because 
the star thinks it is losing negative energy, its density waves get am- 
plified. This intriguing mechanism is similar to that by which spiral 
density waves are produced in galaxies. Although the CFS instability 
was once thought ubiquitous for spinning stars, we now know that neu- 
tron star viscosity will kill it, stabilizing the star and turning off the 
waves, when the star's temperature is above some limit ~ 10 10 K [69] 
and below some limit ~ 10 9 K [70]; and correspondingly, the instability 
should operate only during the first few years of a neutron star's life, 
when 10 9 K £ T £ 10 10 K. 

4. LISA: The Laser Interferometer Space Antenna 

Turn, now, from the high-frequency band, 1 — 10 4 Hz, to the low-frequency 
band, 10~ 4 — 1 Hz. At present, the most sensitive gravitational-wave searches 
at low frequencies are those carried out by researchers at NASA's Jet 
propulsion Laboratory, using microwave-frequency Doppler tracking of in- 
terplanetary spacecraft. These searches are done at rather low cost, piggy- 
back on missions designed for other purposes. Although they have a definite 
possibility of success, the odds are against them. Their best past sensitiv- 
ities to bursts, for example, have been hsB ~ 10 -14 , and prospects are 
good for reaching ~ 10~ 15 — 10~ 16 in the next 5 to 10 years. However, the 
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Figure 7. LISA's orbital configuration. 



strongest low-frequency bursts arriving several times per year might be no 
larger than ~ 10~ 18 ; and the domain of an assured plethora of signals is 
ZiSB ~ 10- 19 -10" 20 . 

In the 2014 time frame, the European Space Agency (ESA) and/or 
NASA is likely to fly a Laser Interferometer Space Antenna (LISA) which 
will achieve hsB & 10~ 20 over the frequency band 3 x 10~ 4 Hz <^ / <5 
3 x l(T 2 Hz. 

4.1. MISSION STATUS 

LISA is largely an outgrowth of 15 years of studies by Peter Bender and 
colleagues at the University of Colorado. Unfortunately, the prospects for 
NASA to fly such a mission have not looked good in the early 1990s. By 
contrast, prospects in Europe have looked much better, so a largely Euro- 
pean consortium was put together in 1993 with Bender's participation but 
under the leadership of Karsten Danzmann (Hannover) and James Hough 
(Glasgow), to propose LISA to the European Space Agency. This proposal 
has met with considerable success; LISA might well achieve approval to fly 
as an ESA Cornerstone Mission around 2014 [71]. Members of the Amer- 
ican gravitation community hope that NASA will join together with ESA 
in this endeavor, and that working jointly, ESA and NASA will be able to 
fly LISA considerably sooner than 2014. 

4.2. MISSION CONFIGURATION 

As presently conceived, LISA will consist of six compact, drag-free space- 
craft (i.e. spacecraft that are shielded from buffeting by solar wind and 
radiation pressure, and that thus move very nearly on geodesies of space- 
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time). All six spacecraft would be launched simultaneously by a single Ar- 
iane rocket. They would be placed into the same heliocentric orbit as the 
Earth occupies, but would follow 20° behind the Earth; cf. Fig. 7. The 
spacecraft would fly in pairs, with each pair at the vertex of an equilateral 
triangle that is inclined at an angle of 60° to the Earth's orbital plane. 
The triangle's arm length would be 5 million km (10 6 times larger than 
LIGO's arms!). The six spacecraft would track each other optically, using 
one- Watt YAG laser beams. Because of diffraction losses over the 5 x 10 6 km 
arm length, it is not feasible to reflect the beams back and forth between 
mirrors as is done with LIGO. Instead, each spacecraft will have its own 
laser; and the lasers will be phase locked to each other, thereby achieving 
the same kind of phase-coherent out-and-back light travel as LIGO achieves 
with mirrors. The six-laser, six-spacecraft configuration thereby functions 
as three, partially independent but partially redundant, gravitational-wave 
interferometers . 

4.3. NOISE AND SENSITIVITY 

Figure 8 depicts the expected noise and sensitivity of LISA in the same 
language as we have used for LIGO (Fig. 3). The curve at the bottom of 
the stippled region is /i rms , the rms noise, in a bandwidth equal to fre- 
quency, for waves with optimum direction and polarization. The top of the 
stippled region is Jisb = ^V^h ims , the sensitivity for high-confidence detec- 
tion (S/N = 5) of a broad-band burst coming from a random direction, 
assuming Gaussian noise. 

At frequencies / P> 10~ 3 Hz, LISA's noise is due to photon counting 
statistics (shot noise). The noise curve steepens at / ~ 3 x 10~ 2 Hz because 
at larger / than that, the waves' period is shorter than the round-trip 
light travel time in one of LISA's arms. Below lCP 3 Hz, the noise is due 
to buffeting-induced random motions of the spacecraft that are not being 
properly removed by the drag-compensation system. Notice that, in terms 
of dimensionless amplitude, LISA's sensitivity is roughly the same as that of 
LIGO's first interferometers (Fig. 3), but at 100,000 times lower frequency. 
Since the waves' energy flux scales as f 2 h 2 , this corresponds to 10 10 better 
energy sensitivity than LIGO. 

4.4. OBSERVATIONAL STRATEGY 

LISA can detect and study, simultaneously, a wide variety of different 
sources scattered over all directions on the sky. The key to distinguish- 
ing the different sources is the different time evolution of their waveforms. 
The key to determining each source's direction, and confirming that it is 
real and not just noise, is the manner in which its waves' amplitude and 
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Figure 8. LISA's projected broad-band noise h m s and sensitivity to bursts hss, com- 
pared with the strengths of the waves from several low-frequency sources. [Note: When 
members of the LISA team plot curves analogous to this, they show the sensitivity curve 
(top of stippled region) in units of the amplitude of a periodic signal that can be detected 
with S/N = 5 in one year of integration; that sensitivity to periodic sources is related to 

the hsB used here by hsp = hsB /y/*3 x 10 7 sec.] 

frequency are modulated by LISA's complicated orbital motion — a motion 
in which the interferometer triangle rotates around its center once per year, 
and the interferometer plane rotates around the normal to the Earth's orbit 
once per year. Most sources will be observed for a year or longer, thereby 
making full use of these modulations. 

5. Low-Frequency Gravitational- Wave Sources 

5.1. WAVES FROM BINARY STARS 

LISA has a large class of guaranteed sources: short-period binary stars in 
our own galaxy. A specific example is the classic binary 44 i Boo (HD133640), 
a 1.35Mq/ O.68M0 system just 12 parsecs from Earth, whose wave fre- 
quency / and characteristic amplitude h c = h^fn are depicted in Fig. 8. 
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(Here h is the waves' actual amplitude and n = fr is the number of wave 
cycles during f =1 year of signal integration.) Since 44 i Boo lies right on 
the hsB curve, its signal to noise ratio in one year of integration should be 
S/N = 5. 

To have an especially short period, a binary must be made of especially 
compact bodies — white dwarfs (WD), neutron stars (NS), and/or black 
holes (BH). WD/ WD binaries are thought to be so numerous that they 
might produce a stochastic background of gravitational waves, at the level 
shown in Fig. 8, that will hide some other interesting waves from view [72]. 
Since WD/ WD binaries are very dim optically, their actual numbers are 
not known for sure; Fig. 8 might be an overestimate. 

Assuming a NS/NS coalescence rate of 1 each 10 5 years in our galaxy 
[6, 7], the shortest period NS/NS binary should have a remaining life of 
about 5 x 10 4 years, corresponding to a gravitational-wave frequency today 
of / ~ 5 x 10~ 3 Hz, an amplitude (at about lOkpc distance) h ~ 4 x 1CT 22 , 
and a characteristic amplitude (with one year of integration time) h c ~ 
2 x 10~ 19 . This is depicted in Fig. 8 at the right edge of the region marked 
"brightest NS/NS binaries" . These brightest NS/NS binaries can be studied 
by LISA with the impressive signal to noise ratios S/N ~ 50 to 500. 

5.2. WAVES FROM THE COALESCENCE OF MASSIVE BLACK HOLES IN 
DISTANT GALAXIES 

LISA would be a powerful instrument for studying massive black holes 
in distant galaxies. Figure 8 shows, as examples, the waves from several 
massive black hole binaries at 3Gpc distance. The waves sweep upward in 
frequency (rightward in the diagram) as the holes spiral together. The black 
dots show the waves' frequency one year before the holes' final collision 
and coalescence, and the arrowed lines show the sweep of frequency and 
characteristic amplitude h c = hy/n during that last year. For simplicity, the 
figure is restricted to binaries with equal-mass black holes: 1O 4 M /1O 4 M , 
1O 5 M /1O 5 M , and 1O 6 M /1O 6 M . 

By extrapolation from these three examples, we see that LISA can study 
much of the last year of inspiral, and the waves from the final collision and 
coalescence, whenever the holes' masses are in the range 3 x 1O 4 M ^ M ^ 
1O 8 M . Moreover, LISA can study the final coalescences with remarkable 
signal to noise ratios: S/N 1000. Since these are much larger S/N's 
than LIGO/VIRGO is likely to achieve, we can expect LISA to refine the 
experimental understanding of black-hole physics, and of highly nonlinear 
vibrations of warped spacetime, which LIGO/VIRGO initiates — provided 
the rate of massive black-hole coalescences is of order one per year in the 
Universe or higher. The rate might well be that high, but it also might be 
much lower. 
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By extrapolating Fig. 8 to lower BH/BH masses, we see that LISA 
can observe the last few years of inspiral, but not the final collisions, of 
binary black holes in the range 100M Q <> M <> 1O 4 M , out to cosmological 
distances. 

Extrapolating the BH/BH curves to lower frequencies using the for- 
mula (time to final coalescence) oc f~ 8 ^ 3 , we see that equal-mass BH/BH 
binaries enter LISA's frequency band roughly 1000 years before their fi- 
nal coalescences, more or less independently of their masses, for the range 
1OOM <, M <, 1O 6 M . Thus, if the coalescence rate were to turn out to 
be one per year, LISA would see roughly 1000 additional massive binaries 
that are slowly spiraling inward, with inspiral rates df /dt readily measur- 
able. From the inspiral rates, the amplitudes of the two polarizations, and 
the waves' harmonic content, LISA can determine each such binary's lumi- 
nosity distance, redshifted chirp mass (1 + z)M c , orbital inclination, and 
eccentricity; and from the waves' modulation by LISA's orbital motion, 
LISA can learn the direction to the binary with an accuracy of some tens 
of arcminutes. 

5.3. WAVES FROM COMPACT BODIES SPIRALING INTO MASSIVE 
BLACK HOLES IN DISTANT GALAXIES 

When a compact body with mass /j, spirals into a much more massive black 
hole with mass M, the body's orbital energy E at fixed frequency / (and 
correspondingly at fixed orbital radius a) scales as E oc fi, the gravitational- 
wave luminosity E scales as E oc fi 2 , and the time to final coalescence thus 
scales as t ~ E/E oc l/fi. This means that the smaller is fi/M, the more 
orbits are spent in the hole's strong-gravity region, a & lOGM/c 2 , and thus 
the more detailed and accurate will be the map of the hole's spacetime 
geometry, which is encoded in the emitted waves. 

For holes observed by LIGO/VIRGO, the most extreme mass ratio that 
we can hope for is fi/M ~ 1M /3OOM , since for M > 3OOM the inspiral 
waves are pushed to frequencies below the LIGO/VIRGO band. This limit 
on fi/M seriously constrains the accuracy with which LIGO/VIRGO can 
hope to map out the spacetime geometries of black holes and test the black- 
hole no-hair theorem (Sec. 3.1). By contrast, LISA can observe the final 
inspiral waves from objects of any mass M ^ O.5M spiraling into holes of 
mass 3 x 1O 5 M £ M £ 3 x 1O 7 M . 

Figure 8 shows the example of a 1OM black hole spiraling into a 1O 6 M 
hole at 3Gpc distance. The inspiral orbit and waves are strongly influenced 
by the hole's spin. Two cases are shown [73]: an inspiraling circular or- 
bit around a non-spinning hole, and a prograde, circular, equatorial or- 
bit around a maximally spinning hole. In each case the dot at the upper 
left end of the arrowed curve is the frequency and characteristic ampli- 
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tude one year before the final coalescence. In the nonspinning case, the 
small hole spends its last year spiraling inward from r ~ 7AGM/c 2 (3.7 
Schwarzschild radii) to its last stable circular orbit at r = 6GM/c 2 (3 
Schwarzschild radii). In the maximal spin case, the last year is spent trav- 
eling from r = 6GM/c 2 (3 Schwarzschild radii) to the last stable orbit 
at r = GM/c 2 (half a Schwarzschild radius). The ~ 10 5 cycles of waves 
during this last year should carry, encoded in themselves, rather accurate 
values for the massive hole's lowest few multipole moments [74]. If the mea- 
sured moments satisfy the "no-hair" theorem (i.e., if they are all determined 
uniquely by the measured mass and spin in the manner of the Kerr metric) , 
then we can be sure the central body is a black hole. If they violate the 
no-hair theorem, then (assuming general relativity is correct), either the 
central body was not a black hole, or an accretion disk or other material 
was perturbing its orbit [75]. From the evolution of the waves one can hope 
to determine which is the case, and to explore the properties of the central 
body and its environment. 

Models of galactic nuclei, where massive holes reside, suggest that in- 
spiraling stars and small holes typically will be in rather eccentric orbits 
[76]. This is because they get injected into such orbits via gravitational 
deflections off other stars, and by the time gravitational radiation reaction 
becomes the dominant orbital driving force, there is not enough inspiral left 
to fully circularize their orbits. Such orbital eccentricity will complicate the 
waveforms and complicate the extraction of information from them. Efforts 
to understand the emitted waveforms are just now getting underway. 

The event rates for inspiral into massive black holes are not at all well 
understood. However, since a significant fraction of all galactic nuclei are 
thought to contain massive holes, and since white dwarfs and neutron stars, 
as well as small black holes, can withstand tidal disruption as they plunge 
toward the massive hole's horizon, and since LISA can see inspiraling bodies 
as small as ~ 0.5M© out to 3Gpc distance, the event rate is likely to be 
interestingly large. 

6. Conclusion 

It is now 35 years since Joseph Weber initiated his pioneering development 
of gravitational-wave detectors [77] and 25 years since Forward and Weiss 
initiated work on interferometric detectors. Since then, hundreds of talented 
experimental physicists have struggled to improve the sensitivities of these 
instruments. At last, success is in sight. If the source estimates described 
in this lecture are approximately correct, then the planned interferometers 
should detect the first waves in 2001 or several years thereafter, thereby 
opening up this rich new window onto the universe. One payoff should 
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be deep new insights into compact astrophysical bodies and their roles in 
binary systems. 
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